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Abstract: In supersymmetric models with gravity-mediated SUSY breaking and gaugino 
mass unification, the predicted relic abundance of neutralinos usually exceeds the strict 
limits imposed by the WMAP collaboration. One way to obtain the correct relic abundance 
is to abandon gaugino mass universality and allow a mixed wino-bino lightest SUSY particle 
(LSP). The enhanced annihilation and scattering cross sections of mixed wino dark matter 
(MWDM) compared to bino dark matter lead to enhanced rates for direct dark matter 
detection, as well as for indirect detection at neutrino telescopes and for detection of dark 
matter annihilation products in the galactic halo. For collider experiments, MWDM leads 
to a reduced but significant mass gap between the lightest neutralinos so that Z2 two-body 
decay modes are usually closed. This means that dilepton mass edges- the starting point 
for cascade decay reconstruction at the CERN LHC- should be accessible over almost all of 
parameter space. Measurement of the — mass gap at LHC plus various sparticle 
masses and cross sections as a function of beam polarization at the International Linear 
Collider (ILC) would pinpoint MWDM as the dominant component of dark matter in the 
universe. 
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1. Introduction 

In supersymmetric models of particle physics, i?-parity is often imposed to avoid too rapid 
proton decay which can be induced by superpotential terms which violate baryon and 
lepton number conservation. One of the byproducts of i?-parity conservation is that the 
lightest supersymmetric particle is absolutely stable, making it a good candidate particle to 
make up the bulk of dark matter (DM) in the universe. In gravity-mediated SUSY breaking 
models, dark matter candidate particles include the lightest neutralino or the gravitino. 
Here we will focus on the lightest neutralino Z\ [jy] ; recent results on TeV scale gravitino dark 
matter can be found in Ref. 0. The relic density of neutralinos in supersymmetric models 
can be straightforwardly calculated by solving the Boltzmann equation for the neutralino 
number density Q. The central part of the calculation is to evaluate the thermally averaged 
neutralino annihilation and co-annihilation cross section times velocity. The computation 
requires evaluating many thousands of Feynman diagrams. Several computer codes are 
now publicly^, |j| available which evaluate the neutralino relic density f^/i 2 . 

The dark matter density of the universe has recently been inferred from the WMAP 
collaboration based on precision fits to anisotropies in the cosmic microwave background 
radiationQ. The WMAP collaboration result for the relic density of cold dark matter 
(CDM) is that 

n CDM h 2 = 0.113 ±0.009. (1.1) 

This result imposes a tight constraint on supersymmetric models which contain a dark 
matter candidate 0. 

Many analyses have been recently performed in the context of the paradigm minimal 
supergravity model[§] (mSUGRA), where the parameter space is given by mo, m^, Aq, tan j3 
and sign(fi). The mSUGRA model assumes the minimal supersymmetric model (MSSM) 
is valid between the mass scales Q = Mqut and Q = Msusy- A common mass mo 
(7714/2) ((A))) is assumed for all scalars (gauginos) ((trilinear soft breaking parameters)) 
at Q = Mgut, while the bilinear soft term B is traded for tan /3, the ratio of Higgs vevs, 
via the requirement of radiative electroweak symmetry breaking (REWSB). REWSB also 
determines the magnitude, but not the sign, of the superpotential Higgs mass term \x. 
Weak scale couplings and soft parameters can be computed via renormalization group 
(RG) evolution from Q = Mqut to Q = M wea k- Once weak scale parameters are known, 
then sparticle masses and mixings may be computed, and the associated relic density of 
neutralinos can be determined. 

In most of mSUGRA parameter space, the relic density turns out to be much 

larger than the WMAP value. Many analyses have found just several allowed regions of 
parameter space: 

• The bulk region occurs at low values of mo and m^H- In this region, neutralino 
annihilation is enhanced by i-channel exchange of relatively light sleptons. The bulk 
region, featured prominently in many early analyses of the relic density, has been 
squeezed from below by the LEP2 bound on the chargino mass m^ > 103.5 GeV, 
and from above by the tight bound from WMAP. 
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• The stau co-annihilation region at low mo for almost any m\i2 value where m^ ~ 
nig , so that f% — Z\ and co-annihilation help to reduce the relic density [flO|] . 

• The hyperbolic branch/focus point (HB/FP) region at large mo ~ several TeV, where 
fi becomes small, and neutralinos efficiently annihilate via their higgsino components|ll] 
This is the case of mixed higgsino dark matter (MHDM). 

• The A-annihilation funnel occurs at large tan j3 values when 2m^ i ~ itia and neu- 



tralinos can efficiently annihilate through the broad A and H Higgs resonances [12]. 



In addition, a less prominent light Higgs h annihilation corridor occurs at low m^flj] and 
a top squark co-annihilation region occurs at particular Aq values when m^ ~ [14|. 

Many analyses have also been performed for gravity-mediated SUSY breaking models 
with non-universal soft terms. Non-universality of SSB scalar masses can 1. pull various 
scalar masses to low values so that "bulk" annihilation via i-channel exchange of light 



scalars can occur [15], or 2. they can bring in new near degeneracies of various sparticles 



with the Z\ so that new co- annihilation regions open up [16, 17, 18], or they can 3. bring 
the value of tha into accord with 2m g so that funnel annihilation can occur [19, 17 1, or 
4. they can pull the value of fi down so that higgsino annihilation can occur [19, 20, 17]. 
It is worthwhile noting that all these general mechanisms for increasing the neutralino 
annihilation rate already occur in the mSUGRA model. Moreover, in all these cases the 
lightest neutralino is either bino-like, or a bino-higgsino mixture. 

If non-universal gaugino masses are allowed, then a qualitatively new possibility arises 
that is not realized in the mSUGRA model: that of mixed wino dark matter (MWDM). 
In this case, if the SU(2) gaugino mass M2 is sufficiently low compared to U(l)y gaugino 
mass Mi, then the Z\ can become increasingly wino-like. The Z\ — W\i — W coupling 
becomes large when Z\ becomes wino-like, resulting in enhanced Z\Z\ — > W + W~ annihi- 
lations. Moreover, coannihilations with the lightest chargino and with the next-to-lightest 
neutralino help to further suppress the LSP thermal relic abundance. 

Non-universal gaugino masses can arise in supersymmetric models in a number of 



ways [21 



In supergravity GUT models, the gauge kinetic function (GKF) /ab must transform 
as the symmetric product of two adjoints. In minimal supergravity, the GKF trans- 
forms as a singlet. In SU(5) SUGRA-GUT models, it can also transform as a 24, 
75 or 200 dimensional representation [p2| , while in 50(10) models it can transform 
as 1, 54, 210 and 770 dimensional representations |23|, E4|]. Each of these non-singlet 
cases leads to unique predictions for the ratios of GUT scale gaugino masses. Fur- 
thermore, if the GKF transforms as a linear combination of these higher dimensional 
representations, then essentially arbitrary gaugino masses are allowed. 

Non-universal gaugino masses are endemic to heterotic superstring models with orb- 



ifold compactification where SUSY breaking is dominated by the moduli fields [25]. 



Additionally, in extra-dimensional SUSY GUT models where SUSY breaking is com- 
municated from the SUSY breaking brane to the visible brane via gaugino mediation, 
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various patterns of GUT scale gaugino masses can occur, including the case of com- 
pletely independent gaugino masses] 26 1. 



In this report, we will adopt a phenomenological approach, and regard the three MSSM 
gaugino masses as independent parameters, with the constraint that the neutralino relic 
density should match the WMAP measured value. 

Much previous work has been done on evaluating the relic density in models with 
gaugino mass non- universality. In AMSB models [27], the Z\ is almost pure wino, so that 
fi^/i 2 as predicted by the Boltzmann equation is typically very low. Moroi and Randall[28] 
proposed moduli decay to wino-like neutralinos in the early universe to account for the 
dark matter density. Already in 1991, Griest and Roszkowski had shown that a wide range 
of relic density values could be obtained by abandoning gaugino mass universality^]. 
Corsetti and Nath investigated dark matter relic density and detection rates in models with 
non-minimal SU(5) GKF and also in O-II string models]^. Birkedal-Hanson and Nelson 
showed that a GUT scale ratio M1/M2 ~ 1.5 would bring the relic density into accord with 
the measured CDM density via MWDM, and also presented direct detection rates pl| . 
Bertin, Nezri and Orloff showed variation of relic density and enhancement in direct and 
indirect DM detection rates as non-universal gaugino masses were varied [32]. Bottino et 
al. performed scans over independent weak scale parameters to show variation in indirect 
DM detection rates, and noted that neutralinos as low as 6 GeV are allowed |33|]. Belanger 
et al. presented relic density plots in the rriQ vs.rri\/2 plane for a variety of universal and 
non-universal gaugino mass scenarios, and showed that large swaths of parameter space 
open up when the SU(3) gaugino mass M3 becomes small [34]. Mambrini and Munoz and 
also Cerdeno and Munoz showed direct and indirect detection rates for model with scalar 
and gaugino mass non-universality [^]. Auto et aZ.pq| used non- universal gaugino masses 
to reconcile the predicted relic density in models with Yukawa coupling unification with the 
WMAP result. Masiero, Profumo and Ullio exhibit the relic density and direct and indirect 
detection rates in split supersymmetry where Mi , and fx are taken as independent weak 
scale parameters with ultra-heavy squarks and sleptons[^6|. 

In this paper, we will adopt a model with GUT scale parameters including universal 
scalar masses, but with independent gaugino masses leading to MWDM. We will assume 
all gaugino masses to be of the same sign. The opposite sign situation leads to a distinct 
DM scenario and will be addressed soon[37[. We will adjust the gaugino masses such that 
Z\ receives just enough of a wino component so that it makes up the entire CDM density 
as determined by WMAP without the need for late-decaying moduli fields. In fact, the 
wino component of the Z\ is usually of order 0.1 — 0.2, so that the Z\ is still mainly 
bino-like, but with a sufficiently large admixture of wino as to match the WMAP result 
on QcDAth 2 . In Sec. |2|, we present the parameter space for MWDM, and show how the 
assumption of MWDM influences the spectrum of sparticle masses. In Sec. ||, we show 
rates for direct and indirect detection of MWDM. These rates are usually enhanced relative 
to mSUGRA due to the enhanced wino component of the Z\. In Sec. [|, we investigate 
consequences of MWDM for the CERN LHC and the international linear e + e~ collider 
(ILC). The goal here is to devise a set of measurements that can differentiate MWDM 
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from the usual case of bino-like DM or MHDM as expected in the mSUGRA model. For 
MWDM, the neutralino mass gap ~ m z 1 ls a l m °st always less than Mz, so that two- 
body decays of Z2 are closed, and three body decays are dominant. The m^ 2 — m z 1 mass 
gap is directly measurable at the CERN LHC via the well-known edge in the m(£ + l~) 
distribution. The correlation of the Z2 — Z\ mass gap against direct and indirect detection 
rates provides a distinction between the possible DM candidates. Measurements of the 
m ~ m z 1 mass § a P a t the LHC combined with measurements of chargino and neutralino 
masses and production cross sections as a function of beam polarization at the ILC would 
provide the ultimate determination of the presence of MWDM in the universe. In Sec. || 
we present our conclusions. 



2. Relic density and sparticle mass spectrum 

Our goal is to explore SUGRA models with non-universal gaugino masses leading to 
MWDM with a neutralino relic density in accord with the WMAP result. To do so, 
we adopt the subprogram Isasugra, which is a part of the Isajet 7.72 event generator 
program|38|]. Isasugra allows supersymmetric spectra generation using a variety of GUT 
scale non-universal soft SUSY breaking terms. The Isasugra spectrum is generated using 
2-loop MSSM RGEs for coupling and soft SUSY breaking term evolution. An iterative ap- 
proach is used to evaluate the supersymmetric spectrum. Electroweak symmetry is broken 
radiatively, so that the magnitude, but not the sign, of the superpotential fi parameter is 
determined. The RG- improved 1-loop effective potential is minimized at an optimized scale 
which accounts for leading 2-loop terms. Full 1-loop radiative corrections are incorporated 
for all sparticle masses. To evaluate the neutralino relic density, we adopt the IsaReD 
program Q, which is based on CompHEP to compute the several thousands of neutralino 
annihilation and co-annihilation Feynman diagrams. Relativistic thermal averaging of the 



cross section times velocity is performed! 39 1. The parameter space we consider is given by 

m , mj/2, A Q , tan/5, sign(n), M x or M 2 , (2.1) 

where we take either M\ or M2 to be free parameters, and in general not equal to m^ 2 . 

In Fig. 1, we show our first result. Here, we take mo = m^2 = 300 GeV, with 
Aq = 0, tan (3 = 10, fj, > with m t = 178 GeV. We plot the neutralino relic density 
fi^/i 2 in frame a) versus variation in the U(l) gaugino mass M\. At M\ = 300 GeV, 
we are in the mSUGRA case, and ^ Zl ^ 2 = ^ ^> so that the model would be excluded 
by WMAP. By decreasing Mi, the bino-like neutralino becomes lighter until two dips 
in the relic density occur. These correspond to the cases where 2m^ ~ and Mz 
as one moves towards decreasing Mi, i.e. one has either light Higgs h or Z resonance 
annihilation. As Mi increases past its mSUGRA value, the Z\ becomes increasing wino- 
like, and the relic density decreases. The W — Wi t 2 — Z\ coupling is proportional to the 
SU(2)l gaugino component of the neutralino, (and also to the Higgsino components), and 
so ^i^i — > W + W~ annihilation becomes enhanced, and the relic density is lowered. In 
this case, the WMAP ^ Zl ^ 2 vame * s reached for Mi = 490 GeV. For still higher Mi values, 
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h 2 drops precipitously, so that other non-neutralino dark matter candidates would have 
to exist to account for the dark matter density in the universe. In frame b), we show the 
bino/wino fraction R^n? of the Z\. Here, we adopt the notation of Ref. ^0|, wherein 
the lightest neutralino is written in terms of its (four component Major ana) Higgsino and 
gaugino components as 

Z x = «f Vfto + 4%° + 4 1)A 3 + v£h , (2.2) 

where R^ = \v% | and = \v± |. While R^ increases as M\ increases, its value when 
Sl^/i 2 reaches the WMAP value is still only ~ 0.25, while ~ 0.9. Thus, the Z\ is still 
mainly bino-like, with just enough admixture of wino to give the correct relic density. This 
corresponds to the case of MWDM. A similar plot is obtained by lowering M2, rather than 
raising Mi, as shown in Fig. 0. 

By raising or lowering the GUT scale gaugino masses in SUGRA models, the mass of 
the neutralinos will obviously change since M\ and M2 enter directly into the neutralino 
mass matrix. However, various other sparticle masses will also be affected by varying the 
gaugino masses, since these feed into the soft term evolution via the RGEs. In Fig. [3L 
we show the variation of the sparticle mass spectrum with respect to the GUT scale ratio 
Mi/mi/2 for the same parameters as in Fig. |]. When Mi/rrii/2 = 1, there is a relatively 
large mass gap between Z2 and Z\\ m^ 2 — = 106.7 GeV. As Mi is increased until 
fi^/i 2 = 0.11, the mass gap shrinks to fn^ 2 — m z 1 = 31-9 GeV. The light chargino mass 
remains essentially constant in this case, since M2 remains fixed at 300 GeV. However, 
we notice that as Mi increases, the e^j, jlji and fi masses also increase, since Mi feeds into 
their mass evolution via RGEs. As the coefficient appearing in front of Mi in the RGEs is 
larger (and with the same sign) for the right handed sfermions than for the left handed ones, 
one expects, in general, a departure from the usual mSUGRA situation where the lightest 
sleptons are right-handed. As a matter of fact, whereas in mSUGRA mg L >> rrie R , in the 
case of MWDM, instead, for the particular parameter space slice under consideration, we 




Figure 1: A plot of a) relic density flcDMh 2 and b) bino/wino component of the lightest neutralino 
as a function of M\ for mo = 300 GeV, mi/2 = 300 GeV, Aq = 0, tan f3 = 10, fi > and mt = 178 
GeV. 
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find that ~ rrie R - As shown in the figure, the right-handed squark masses also increase 
with increasing Mi, although the relative effect is less dramatic than the case involving 
sleptons: the dominant driving term in the RGEs is, in this case, given by M3 (absent in 
the case of sleptons), hence variations in the GUT value of Mi produce milder effects. 

In Fig. |4|, we show a plot of sparticle masses for the same parameters as in Fig. [|, but 
versus Mijrai^- In this M2 is decreased from its mSUGRA value of 300 GeV, 

the W\ and Z2 masses decrease until £1% h 2 reaches 0.11, where now m^ 2 — = 22.9 
GeV. In this case, with decreasing M2, the left- slepton and sneutrino masses also decrease, 
again leading to ~ rrie R . The left-handed squark masses similarly decrease. Right- 
handed sfermion masses are, instead, not affected, with the net result that the mSUGRA 
rrig L » m e - fi hierarchy is again altered. 

The effect of varying gaugino masses on the allowed region of parameter space is 
illustrated in Fig. [|. Here, in frame a), we show the case of the mSUGRA model in 
the mo vs. mi/2 plane for Aq = 0, tan/? = 10 and \i > 0. The red shaded regions are 
disallowed by either a stau LSP (left side of plot) or lack of REWSB (lower edge of plot). 
The blue shaded region has a chargino with mass m^ < 103.5 GeV, thus violating bounds 
from LEP2. The dark green shaded region has 0.094 < ^^ ± h 2 < 0.129, in accord with the 
WMAP measurement. The light green shaded region has f^fr- 2 < 0.094, so that additional 
sources of dark matter would be needed. We see the stau co-annihilation region appearing 
along the left edge of the allowed parameter space, and the bulk region appearing at low mo 
and low mi/2- The h annihilation corridor appears also at low mi/2 along the edge of the 
LEP2 excluded region. In frame b), we take Mi/mi/2 = 1.5, so that the Z\ becomes more 
wino-like. In response, we see that a large new bulk region has appeared at low mo and 
low mx/2- In frame c), we increase M\/mi/2 to 1.75. In this case, most of the mo vs. my 2 
plane is now allowed, although much of it has ^^Ji 2 below the WMAP central value for 
the CDM relic density. 



m =300 GeV, m 1/2 =300 GeV, tan(3=10, A Q =0, u>0, m =178 GeV 




M 2 (GeV) M 2 (GeV) 



Figure 2: A plot of a) relic density flcDMh 2 and 6) bino/wino component of the lightest neutralino 
as a function of Mi for mo = 300 GeV, m\ii — 300 GeV, Aq = 0, tan j3 = 10, fi > and mt = 178 
GeV. 
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It should be apparent now that any point in the rriQ vs. m\/2 plane can become WMAP 
allowed by either increasing M\ or decreasing M2 to a suitable degree as to obtain MWDM. 
To illustrate this, we plot in Fig. |6] the ratio T\ = Mi/m^ in frame a) or T2 = Milvnyji 
in frame b) needed to achieve a relic density in accord with the WMAP central value. We 
see in frame a) that r\ increases as one moves from lower-left to upper-right, reflecting the 
greater wino component of Z\ that is needed to overcome the increasing Q^^ 2 which is 
expected in the mSUGRA model. We also see on the left side of the plot that r\ < 1 is 
allowed, since then Q^h 2 < 0.11 already in the mSUGRA case. The structure at high m-y^ 
and mo ~ 400 — 500 GeV results because increasing M\ increases until 2m^ ~ tua 
and the A- funnel begins to come into effect (even though tan/3 is small). 

3. Direct and indirect detection of mixed wino dark matter 

In this section, we turn to consequences of MWDM for direct and indirect detection of 
neutralino dark matter ^l], E^]. We adopt the DarkSUSY code[^3|], interfaced to Isajet, for 
the computation of the various rates, and resort to the Adiabatically Contracted N03 Halo 
model p4}| for the dark matter distribution in the Milky Way 1 . We evaluate the following 
neutralino DM detection rates: 



Direct neutralino detection via underground cryogenic detectors [48]. Here, we com- 
pute the spin independent neutralino-proton scattering cross section, and compare 
it to expected sensitivities @ for Stage 2 detectors (CDMS2@, Edelweiss2[|l|, 
CRESST2||, ZEPLIN2||) and for Stage 3, ton-size detectors (XENON@, Genius|5 
ZEPLIN4[^] and WARP [57|1 ) . We take here as benchmark experimental reaches of 
Stage 2 and Stage 3 detectors the projected sensitivities of, respectively, CDMS2 and 
XENON 1-ton at the corresponding neutralino mass. 

Indirect detection of neutralinos via neutralino annihilation to neutrinos in the core of 



the Sun [ 58 1 . Here, we present rates for detection of —* [i conversions at Antares[59] 
or IceCube|30|]. The reference experimental sensitivity we use is that of IceCube, with 
a muon energy threshold of 25 GeV, corresponding to a flux of about 40 muons per 
km 2 per year. 

Indirect detection of neutralinos via neutralino annihilations in the galactic center 
leading to gamma rays|3l]], as searched for by EGRET |62|, and in the future by 



GLAST[63|. We evaluate the integrated continuum 7 ray flux above a E 7 = 1 GeV 
threshold, and assume a GLAST sensitivity of l.OxlO -10 cm~ 2 s _1 . 

Indirect detection of neutralinos via neutralino annihilations in the galactic halo 
leading to cosmic antiparticles, including positrons |64| (HEAT f65j| , Pamela []66| and 
AMS-02|67|]), antiprotons@ (BESS§9|, Pamela, AMS-02) and anti-deuterons (Ds) 
(BESS|7(|, AMS-02, GAPS[71|). For positrons and antiprotons we evaluate the av- 



eraged differential antiparticles flux in a projected energy bin centered at a kinetic 



For a comparison of the implications of different halo model choices for indirect DM detection rates, 
see e.g. Refs. || || 0, 0. 
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energy of 20 GeV, where we expect an optimal statistics and signal-to-background ra- 
tio at space-borne antiparticles detectors [47, ff2^ . We use as benchmark experimental 
sensitivity that of the Pamela experiment after three years of data-taking. Finally, 
the average differential antideuteron flux has been computed in the 0.1 < < 0.4 
GeV range, and compared to the estimated GAPS sensitivity[71]. 



In Fig. [7], we show various direct and indirect DM detection rates for tuq = = 300 
Gev, with Aq = 0, tan/3 = 10 and /i > 0, while M\ is allowed to vary. The M\ value 
corresponding to the mSUGRA model is denoted by a dot-dashed vertical line, while the 
one where Q^h 2 = 0.11 by a dashed vertical line denoted WMAP. 

In frame a), we plot the spin-independent neutralino-proton scattering cross section. 
Both the squark-mediated and Higgs mediated neutralino-proton scattering amplitudes are 
enhanced by more than one order of magnitude due to the increasing wino nature of the Z\ . 
The reason for the enhancement is traced back to the structure of the neutralino-quark- 
squark and neutralino-neutralino-Higgs couplings, where the wino fraction is weighed by 
the SU(2) coupling, while the bino fraction by the (smaller) U(l) coupling. 

In frame b), we show the flux of muons from neutralino pair annihilations in the core 
of the Sun. While the muon flux is below the reach of IceCube in the mSUGRA case, it 
has climbed into the observable region when the Z\ has become sufficiently wino-like as to 
fulfill the WMAP measured DM relic density. 

In frames c), d), e) and f) we show the flux of photons, positrons, antiprotons and 
antideuterons, respectively. The results here are plotted as ratios of fluxes normalized to the 
mSUGRA point, in order to give results that are approximately halo-model independent. 
(We do show the above described expected experimental reach lines as obtained by using 
the Adiabatically Contracted N03 Halo model[^4|].) All rates are enhanced, with respect 
to the mSUGRA case, by 2 to 3 orders of magnitude, due to the increasing cross section 
for Z\Z\ — > W + W~ annihilation in the galactic halo. In particular, antimatter fluxes are 
always below future sensitivities for the mSUGRA setup, while they all rise to a detectable 
level when the WMAP point is reached. 

In Fig. we show the same direct and indirect DM detection rates as in Fig. [7|, 
except this time versus M<i instead of M\ . In this case, the various rates are all increasing 
as Mi decreases, entering the region of MWDM. Indirect detection rates again feature 
enhancements as large as 2 orders of magnitude with respect to the mSUGRA scenario, 
when the WMAP relic abundance is reached. The abrupt decrease in the rates below 
Mi ~ 100 GeV is due, instead, to the < mw threshold. 

In Fig. ^, we show regions of the mo vs. mi/2 plane for Aq = 0, tan/3 = 10 and 
fi > which are accessible to various direct and indirect DM search experiments. The 
visibility criteria we adopt here follow the same approach outlined in Ref. [17|. The gray 



shaded regions in the plots are already excluded, at 95% C.L., by a x 2 analysis of the 
computed signal plus background p flux compared to the available antiprotons data (for 
details see[fl7|). Observable rates for 7 detection by GLAST occur throughout all three 
planes, due to the high DM density assumed at the galactic core in the N03 halo model. 
In frame a), we show the case of the mSUGRA model. Only small regions at low uiq 
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and low m-i/2 are accessible to D searches by GAPS and p searches by Pamela. A tiny 
region is also accessible to CDMS2, and a much larger region is accessible to Stage 3 direct 
detection experiments such as XENON. In frame b), we increase Mi(Mgut) at every point 
in the plane as in Fig. |6| until Sl^^- 2 = 0.11. The corresponding neutralino masses are 
therefore accordingly increased with respect to the mSUGRA case. Nevertheless, we see 
that the regions accessible to direct and indirect DM detection have vastly increased. The 
D search by GAPS can cover < 400 — 500 GeV. The e + and p searches by Pamela 

can see to m 1 / 2 ~ 250 GeV and 350 GeV, respectively. In addition, a region has opened 
up which is accessible to IceCube searches for dark matter annihilation in the core of the 
Sun. The Stage 3 dark matter detectors can see most of the tjiq vs. m\/ 2 plane, with the 
exception of the region at large m^j and low tuq where a much lower wino component of 
the Z\ is required to bring the relic density into line with the WMAP measurement (here, 
early universe Z\ Z\ annihilations are already somewhat enhanced by the proximity of the 
vl-pole and the stau co-annihilation region). In frame c), we show again the mo vs. m\/ 2 
plane, but this time we have reduced M 2 until the Q^^ 2 = 0.11 value is reached. Again, 
many of the direct and indirect detection regions are expanded compared to the mSUGRA 
case. We remark that, although in this last case the neutralino mass is lower than in the 
case shown in frame b), direct detection and neutrino fluxes are somewhat less favored. 
This depends on the relative higgsino fraction, which critically enters in the neutralino- 
proton scattering cross section as well as in the neutralino capture rate in the Sun: raising 
Mi shifts the gaugino masses closer to (i, hence increasing the higgsino fraction and the 
resulting neutralino cross sections off matter. 

4. Mixed wino dark matter at colliders 

An important question is whether collider experiments would be able to distinguish the 
case of MWDM from other forms of neutralino DM such as bino-DM or MHDM as occur in 
the mSUGRA model. We have seen from the plots of sparticle mass spectra that the squark 
and gluino masses vary only slightly with changing M\ or M 2 . However, the chargino and 
neutralino masses change quite a bit, and in fact rather small mass gaps ~ m z l anc ^ 
m_g — nig are in general expected in the case of MWDM as compared with the case from 
models containing gaugino mass unification. 

In Fig. pi], we show contours of the mass gap m^ 2 — mg in the tuq vs. m 1 / 2 plane for 
A = 0, tan (3 = 10 and fj, > for a) the mSUGRA model, b) the case of MWDM where 
Ml is raised at every point until Qg^ 2 —►0.11 and c) the case of MWDM where M 2 is 
lowered until Q,g h 2 — ► 0.11. In the case of the mSUGRA model, most of the parameter 
space has — rn^^ > 90 GeV, which means that Z 2 — ► Z±Z° decay is allowed. When 
this decay is allowed, its branching fraction is always large, unless it competes with other 
two-body decays such as Z 2 — » Z\h or Z 2 — > ff or // (where / is a SM fermion). In the 
case of MWDM in frames b) and c), we see that (aside from the left-most portion of frame 
b), which is not a region of MWDM), the mass gap is much smaller, so that two-body 
decays of Z 2 and W\ are closed and three-body decays are dominant. 



-9- 



When the decays Z2 — ► tl, li — ► Z\lt or Z2 — > Z\U are open (I = e or yu), then 
prospects are good for measuring the Z2 — -Zi mass gap at the CERN LHC and possibly 
at the Fermilab Tevatron. If Z2's are produced at large rates either directly or via gluino 
or squark cascade decays^], it should be possible to identify opposite sign/ same flavor 
dilepton pairs, to reconstruct their invariant mass, and extract the upper edge of the 
invariant mass distribution ]74|]. In Fig. [ll], we show the branching fraction BF(Z2 — > 
Z\e + e~) versus M\ (left-side) or versus M2 (right-side) for a variety of choices of mo, mi/2 
and tan/3. The mSUGRA model value is denoted by the dot-dashed vertical line, while 
the M%,2 value at which ^l^Ji 2 — > 0.11 is indicated by the dotted vertical line. As one 
moves to higher M\ (or lower M2) values, in most cases the leptonic three-body decays of 
Z2 become enhanced, usually because as Mi grows (M2 decreases), the two-body decay 
modes become kinematically closed, and only three-body decays are allowed. Thus, while 
the mSUGRA model yields large rates for Z2 — ► Z\e + e~ only when ~ 220 GeV, this 
decay mode is almost always open in the case of MWDM. The only exception occurs when 
the stau co-annihilation or the >l-funnel act to lower the relic density, so that a large M\ 
or small M2 is not needed to obtain the correct relic density; this, however, is not the case 
of MWDM. 

4.1 CERN LHC 

If the i?-parity conserving MSSM is a good description of nature at the weak scale, then 
multi-jet plus multi-lepton plus Et events should occur at large rates at the CERN LHC, 
provided that m g ~ 2 — 3 TeV. The LHC reach for SUSY in the mSUGRA model has 
been calculated in Ref. IfrEfl . The mSUGRA reach results should also apply qualitatively 
to the MWDM case, since the values of m g and nig change little in going from mSUGRA 
to MWDM, and the reach plots mainly depend on these masses. 

For SUSY searches at the CERN LHC, Hinchliffe et al. have pointed out [[?(]] that an 
approximate value of rriq or rrig can be gained by extracting the maximum in the M e ff 
distribution, where M eff = E T + E T (jet 1) + E T (jet 2) + E T (jet 3) + E T (jet 4). This 
statement holds true in models with MWDM, as well as in models with gaugino mass 
unification, so that the approximate mass scale of strongly interacting sparticles will be 
known soon after a supersymmetry signal has been established. 

In mSUGRA, a dilepton mass edge should be visible in SUSY signal events only if 
m i/2 ~ 250 GeV or if Z2 — ► U, it decays are allowed. In the case of MWDM, the dilepton 
mass edge should be visible over almost all parameter space. We illustrate the situation for 
four case studies listed in Table ||. The first case, labeled mSUGRA, has mo = mjft = 300 
GeV, with A§ = 0, tan/3 = 10 and \x > 0. In this case, gg, gq and qq production occurs 
with a combined cross section of about 12 pb, while the total SUSY cross section is around 
13.4 pb (the additional 1.4 pb comes mainly from -ino pair production and -ino-squark 
or -ino-gluino associated production). The case of MWDM1, with M\ = 490 GeV, has 
similar rates of sparticle pair production. The case of MWDM2, with lighter chargino 
and neutralino masses, has a total production cross section of 19.2 pb, wherein strongly 
interacting sparticles are pair produced at similar rates as in mSUGRA or MWDM1, but 
-ino pairs are produced at a much larger rate ~ 6.1 pb. We also show a case of MHDM 
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parameter 


mSUGRA 


MWDM1 


MWDM2 


MHDM 


Mi 


300 


490 


300 


300 


M 2 


300 


300 


187 


300 


M 


409.2 


410.1 


417.8 


166.1 


m~g 


732.9 


732.8 


733.0 


854.6 


mu L 


720.9 


721.1 


706.9 


3467.2 


m h 


523.4 


526.0 


533.2 


2075.8 


m ~ bl 


650.0 


648.9 


640.2 


2847.0 


m- eL 


364.7 


371.7 


330.0 


O A A f\ T 

3449.7 


rrie R 


QOO Q 

ozz.o 


OOO. ( 


goo *7 
OZZ. t 


OA ACi A 


m W 2 


/I Q 

4oz.y 


A VX H 
400.0 


A "3^ Q 


988 Q 

Zoo.y 


m m 


ooq n 


oo/i n 

ZZ4.U 


loo.o 


i Aa a 
140. 


m z 4 


433.7 


/lor T 

435.7 


436. z 


296.9 


m z 3 


414.8 


415.6 


424.1 


179.0 


m z 2 


223.7 


225.4 


138.8 


159.2 




117.0 


193.5 


115.9 


101.5 


m A 


538.6 


544.1 


523.6 


3409.9 


m H + 


548.0 


553.5 


533.1 


3433.3 


m h 


115.7 


115.8 


115.3 


118.9 




1.3 


0.11 


0.11 


0.13 


BF(b -► si) 


3.2 x 10~ 4 


3.2 x 10~ 4 


3.3 x 10~ 4 


3.4 x 10~ 4 


Aa M 


12.1 x 10~ 10 


11.8 x 10~ 10 


15.9 x 10~ 10 


3.9 x 10~ n 


a sc (Zip) 


2.6 x 10 -8 pb 


2.2 x 10~ 7 pb 


7.1 x 10" 8 pb 


1.8 x 10" 8 pb 



Table 1: Masses and parameters in GeV units for mSUGRA, MWDM and MHDM models. In 
the first three cases, mo = mi/2 = 300 GeV, Aq = 0, tan/3 = 10 and m* = 178 GeV. The case of 
MHDM has the same parameters, except mo = 3451.8 GeV, with m t — 175 GeV. 



from the HB/FP region of the mSUGRA model as an alternative low Z2 — Z\ mass gap 
model to compare against MWDM scenarios. 

We have generated 50K LHC SUSY events for each of these cases using Isajet 7.72, and 
passed them through a toy detector simulation. The toy detector is divided into calorimeter 
cells of size A77 x A(f> = 0.05 x 0.05 extending out to | ?y| < 5, with no transverse shower 
spreading. We invoke EM smearing with 3%/\^E+ .5%, hadronic smearing with 80%/v / £'+ 
3% out to \rj\ = 2.6, and forward calorimeter hadronic smearing with 100%/\/^ + 5%. Jets 
are clustered using a UAl type algorithm with cone size R = \J Arf + A<p 2 = 0.7, with 
Ej e t(min) = 25 GeV. Leptons {I = e or n) with En > 10 GeV are classified as isolated if 
Ex(cone) < 5 GeV in a cone of R = 0.3 about the lepton's direction. Since gluino and 
squark masses of the three case studies are similar to those of LHC point 5 of the study of 
Hinchliffe et aZ.j76|, we adopt the same overall signal selection cuts which gave rise to only 
a small background contamination of mostly signal events: Et > max(100 GeV, 0.2M e ff), 
at least four jets with Et > 50 GeV, where the hardest jet has Et > 100 GeV, transverse 
sphericity S T > 0.2 and M eff > 800 GeV. 
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In these events, we require at least two isolated leptons, and then plot the invariant 
mass of all same flavor /opposite sign dileptons. The results are shown in Fig. [l^. In 
the case of the mSUGRA model, frame a), there is a sharp peak at m(£ + £~) ~ Mz, 
which comes from Z2 — > Z1Z decays where Z2 is produced in the gluino and squark 
cascade decays. In the case of MWDM1 in frame b), we again see a Z° peak, although 
here the Z°s arise from Z3, Z^ and W2 decays. We also see the continuum distribution 
in m(l + l~~) < m^ 2 — = 31.9 GeV. The cross section plotted here is ~ 0.05 pb, 
which would correspond to 5K events in 100 fb _1 of integrated luminosity (the sample 
shown in the figure contains just 406 events). In frame c)- with a cross section of ~ 0.05 
pb (but just 267 actual entries)- we see again the Z° peak, but also we see again the 
m(l + l~) < m z 2 ~ m Zi = 22-9 GeV continuum. In both these MWDM cases, the m z 2 ~ m Zi 
mass edge should be easily measurable. It should also be obvious that it is inconsistent 
with models based on gaugino mass unification, in that the projected ratios Mi : M2 : M3 
will not be in the order 1 :~ 2 :~ 7 as in mSUGRA. Although the Z2 — Z\ mass edge will 
be directly measurable, the absolute neutralino and chargino masses will be difficult if not 
impossible to extract at the LHC. 

In frame d), we show the spectrum from MHDM in the HB/FP region of the mSUGRA 
model. In this Z2 — Z\ mass edge at 57.7 GeV should be visible. It will be 

accompanied by other continuum contributions, since in the case of MHDM with a small 
\i parameter, the Z3, Z4 and W2 should all be relatively light as well. 

4.2 Linear e + e collider 

At a y/s = 500 GeV ILC, the new physics reactions for the four case studies shown in Table 
[l] would include Zh, W^~W^~ , Z1Z2 and Z2Z2 production. It was shown in Ref. |77j that, in 
the case of a small W\—Z\ mass gap, chargino pair production events could still be identified 
above SM backgrounds. The chargino and neutralino masses can be inferred from the 



resultant dijet distribution in W^W^ — * (£ugZi) + (qqZi) events[78, 79. 77|. Alternatively, 
the chargino mass may be extracted from threshold cross section measurements when the 
CM energy of the accelerator is tuned to operate just above e + e~ — ► W^W^ threshold. 
These measurements should allow the absolute mass scale of the sparticles to be pinned 
down, and will complement the Z2 — Z\ mass gap measurement from the CERN LHC. The 
combination of , , and ~ m z l measurements will point to whether or not 
gaugino mass unification is realized in nature. 

In addition, the W^~W^~ , Z1Z2 and Z2Z2 production cross sections can all be measured 
as a function of beam polarization at the ILC. In the mSUGRA model, since W\ and Z2 
are mainly wino-like, they will be produced at high rates for left-polarized electron beams, 
but at low rates for right-polarized beams|79|. The Z1Z2 production cross section also 



has a significant rise to it as beam polarization parameter Pi(e~) is increased from -1 to 
+1. These cross sections are plotted in frame a) of Fig. [13]. In frame b), we show the 
same cross sections, except this time for the case of MWDM1. The W\ is still mainly 
wino-like, and so has a steeply rising cross section as Pi,(e~) is increased. However, in this 
case Z\ and Z2 both have non-negligible bino components, which enhances their couplings 
to right-polarized electrons. Thus, a(e + e~ — > Z\Z<i) in the case of MWDM is a falling 



- 12 - 



distribution vs. Pz,(e~). This is in fact borne out in frame b), and would be a strong signal 
for MWDM! In frame c), we plot the corresponding cross sections for the case of MWDM2. 
Again, Z1Z2 has a (slightly) falling cross section versus Px,(e~), indicating once again the 
presence of MWDM. In frame d), we show the corresponding cross sections for the case of 
MHDM. In this case, numerous other reactions such as W^~W 2 ~ , Z\Z% and Z 2 Z^ should 
likely be kinematically accessible, and their presence will help serve to distinguish MHDM 
from MWDM. 

While a combination of mass measurements at LHC and ILC would help to pin down 
the properties of MWDM, it is worth considering whether the case of MWDM can be 
confused with the case of MHDM, such as occurs in the HB/FP region of the mSUGRA 
model, or in models with non-universal Higgs masses [19, 17]. To answer this, we plot in 
Fig. 14 the Z2 — Z\ mas gap versus for MWDM scenarios which yield Sl^^i 2 = 0.11, 
against mSUGRA models in the HB/FP region which also give Q^^ 1 = 0.11. We see 
that the MWDM points can span the entire range of values shown, but that their 
Z<i — Z\ mass gap is generally of order 15-40 GeV. Models with higher mass gaps are 
usually due to an overlap of MWDM with stau co-annihilation or ^4-funnel annihilation. 
The general trend for — in the MWDM scenario is dictated by the interplay of 
wino coannihilations and of the growing wino component, both functions of the mass gap, 
which suppress to the required level. In contrast, the Z2 — Z\ mass gap associated 

with MHDM in the HB/FP region is generally or order 40-80 GeV, at least until very 
large values of ~ 600 GeV are generated. The largest mass gaps appear beyond the 
top quark mass threshold, whose effect is greatly enhanced, with respect to the MWDM 
case, due to Z and Higgs s-channel exchanges. At larger neutralino masses, the Z2 — Z\ 
mass gap for MHDM shrinks to lower values, since a larger and larger higgsino fraction 
and stronger neutralino/ char gino coannihilations are needed to fulfill the WMAP bound. 
Eventually, a pure higgsino LSP (with — of the orders of few GeV) is needed to 
give tt^h 2 = 0.11, at ~ 1 TeV. For ~ 600 - 800 GeV, the MWDM and MHDM 
Z2 — Z\ mass gaps overlap. In the large mass case, however, the two scenarios could still 
be differentiated by the remaining sparticle mass spectrum (e.g. Z3 would be light in the 
case of MHDM and heavy in the case of MWDM) and by the dependences of cross sections 
on electron beam polarization (if an energetic enough e + e~ collider is built!). 



5. Conclusions 

In this paper, we have considered the phenomenological consequences of mixed wino dark 
matter. MWDM occurs in models with gaugino mass non-universality. MWDM may be 
obtained by modifying the paradigm mSUGRA model by either increasing the GUT scale 
value of M\ or by decreasing M2 until a sufficiently wino-like LSP is obtained as to fulfill 
the WMAP measured value of QcDMh 2 ~ 0.11. If DM in nature is indeed composed of 
MWDM, then a number of consequences occur. In the sparticle mass spectrum, the Z2 — Z\ 
and W\ — Z\ mass gaps are expected to be reduced compared to what is expected in models 
with gaugino mass unification and a large [i parameter. Also, left- and right- sleptons are 
expected to be more nearly mass degenerate. 
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If MWDM comprises the dark matter of the universe, then both direct and indirect 
dark matter detection rates are expected to be enhanced compared to expectations from 
the mSUGRA model. However, to really pinpoint the existence of a partially wino-like 
Z\, collider experiments will be needed. The CERN LHC should be able to measure 
approximately the value of mg, and in MWDM scenarios, also the Z2 — Z\ mass gap from 
the dilepton spectrum from Z^ — > IIZ\ decay. These measurements should be enough to 
establish whether gaugino mass unification holds. Ultimately, a linear e + e~ collider, the 
ILC, operating above W^W^~ and Z1Z2 thresholds will be needed. The ILC should be 
able to measure the absolute W\, Z\ and Z^ masses. The dependence of the associated 
production cross sections on the electron beam polarization will point conclusively to the 
existence of MWDM. 
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Figure 5: WMAP allowed regions in the too vs. TO;iy 2 plane for tan f3 = 10, Aq — 0, [i > and a) 
M 1 /m 1/2 = 1, 6) Mi/to 1/2 = 1.5 and c) Mi/to 1/2 = 1.75. 




Figure 6: Contours of a) ri and b) r 2 in the too vs. mi/ 2 plane for tan/3 = 10, A = 0, /i > 0. 
Each point has f^/i = 0.11. 
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Figure 7: Rates for direct and indirect detection of neutralino dark matter vs. Mi for mo = 
m i/2 — 300 GeV, with tan/3 = 10, Aq = 0, [i > 0. Frames c) -f) show the ratio of indirect detection 
rates compared to the mSUGRA model. In this plot, we adopt the N03 distribution for halo dark 
matter. 
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Figure 8: Rates for direct and indirect detection of neutralino dark matter vs. Mi for mo = 
m 1 / 2 = 300 GeV, with tan/3 = 10, A = 0, \i > 0. Frames c) -f) show the ratio of indirect detection 
rates compared to the mSUGRA model. In this plot, we adopt the N03 distribution for halo dark 
matter. 
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Figure 9: Regions of visibility for direct and indirect dark matter searches in the mo vs. mi/ 2 
plane for tan/3 = 10, Aq — 0, /i > 0. The upper frame a) shows the mSUGRA model, while frame 
b) corresponds to the MWDM model with non- universal Mi and frame c) with non-universal M^. 
In this plot, we adopt the Adiabatically Contracted N03 Halo Model for the galactic dark matter 
distribution. For this halo model, detection of 7s by GLAST should occur over all three planes. 
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mSUGRA: tan$=10, A =0, n >0, m t =178 GeV 
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Figure 10: Contours of — mass gap in the tuq vs. m 1 / 2 plane for tan/3 = 10, Aq = 0, 
fi > and a) mSUGRA model, b) My > m 1/2 MWDM and c) M 2 < m 1/2 MWDM. 
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Figure 11: The branching fraction for — > Z\e + er decay is plotted vs. M\ (left-side) or 
M2 (right-side) for various points in the MWDM model parameter space. The value from 

mSUGRA is denoted by the dot-dashed lines, while the Mi,2 value which gives f^/i 2 = 0.11 is 
indicated by dotted lines. 



- 26 - 



-2 

x 10 



0.007 
$ 0.006 
| 0.005 
"° 0.004 
0.003 
0.002 
0.001 


x 10 



I 



I I I I I I I I I I I I I I 


i 


o.r c 






0.09 I 






0.08 I 






0.07 I 


[_ a) mSUGRA 




0.06 j 






0.05 I 






0.04 I 






0.03 I 






0.02 I 






0.07 I 


i i-nir-J-Hi jn-m-r-LT^ Fl i-h-H-. — hJ p-t-i i-t— 




E 





-2 



0.7 
0.09 
0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 




20 



40 



60 



C) MWDM2 



If 



I I I I I L_U I lL 



n 



80 100 



x 10 





-2 



1\J 



b) MWDM1 



n 



\ 



_J I I I T I | | | | | L 



20 40 60 80 100 



20 40 60 80 100 
m(l + l)(GeV) 




60 80 100 

m(l + l)(GeV) 



Figure 12: Distribution of same flavor/opposite sign dileptons from SUSY events at the CERN 
LHC from a) mSUGRA, b) MWDM1, c) MWDM2 and d) MHDM cases as in Table 0. 
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Figure 13: Plot of cross section for e + e~ — ► W^~W[~ , Z\Zi and Z2Z1 versus electron beam 
polarization Pt,(e") for a = 500 GeV ILC for a) mSUGRA, b) MWDM1, c) MWDM2 and d) 
MHDM with parameters as in table |l|. 
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Figure 14: Correlation between — ra^ and in models with MWDM and MHDM in the 
HB/FP region. 
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